Mean-field-like behavior of the generalized voter-model-class kinetic Ising model 
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We analyze a kinetic Ising model with suppressed bulk noise which is a prominent representative 
of the generalized voter model phase transition. On the one hand we discuss the model in the 
context of social systems, and opinion formation in the presence of a tunable social temperature. 
On the other hand we characterize the abrupt phase transition. The system shows non-equilibrium 
dynamics in the presence of absorbing states. We slightly change the system to get a stationary 
state model variant exhibiting the same kind of phase transition. Using a Fokker-Planck description 
and comparing to mean field calculations, we investigate the phase transition, finite size effects and 
the effect of the absorbing states resulting in a dynamic slowing down. 



I. INTRODUCTION 

Opinion formation models sparked considerable inter- 
est in the physics community, partly due to its close rela- 
tionship to spin models (see [1] for an overview). Simple 
mathematical rules for the outcome of discussions among 
agents determine the dynamics of a system of agents, as 
for example agreement including the time needed for con- 
sensus or disagreement, as well as the spatial spreading 
of opinions in terms of coarsening or segregation. Three 
interesting lines of research in this field are: the voter 
model, the universality class of the generalized voter 
model, and the Sznajd model. 

The Voter Model (VM) [2] on regular lattices describes 
agents with two possible opinions, denoted as spin values 
dzl. The agents are randomly chosen to adopt the opinion 
of one of their nearest neighbors. This parameter free 
Z2 symmetric model includes absorbing states with total 
agreement and thus we have a non-equilibrium system. 
The VM turned out to be one of the rare analytically 
solvable non-equilibrium models (see [1] and references 
therein). Additional interesting properties of the model 
are the lack of surface tension and a slow domain growth 
with diffusively roughened interfaces in two dimensions 

i. 

The universality class of the generalized voter model 
(GVM) [3] is characteristic for systems with parameter- 
ized interactions which show a special non-equilibrium 
phase transition between order (consensus) and disorder 
for certain parameter variations. These systems may in- 
clude the VM at the critical point, and exhibit an abrupt 
phase transition (with a jump in the order parameter 
at the critical parameters), however, show critical di- 
vergences with specific critical exponents for susceptibil- 
ity (7 = 1) and correlation length {v = 1/2) [3]. Be- 
sides the directed percolation phase transition the GVM 
phase transition is a central universality class of non- 
equilibrium phase transitions 0, (H. Many new models 
showing a generalized voter transition were defined and 
investigated, for example using backward Fokker-Planck 
equations, mean field calculations, and Langevin descrip- 
tion [iil. 

The Sznajd model finally puts more emphasis on the 
social interpretation, where persuasiveness increases with 



the number of proponents [9| (see f]\ for an overview of 
dynamics and model variants). A similar effect of win- 
ninglocal majority is incorporated in Majority Rule mod- 
els 111, [loj . A model with agents following local as well as 
global majorities is used in the context of stock markets 



We here analyze a kinetic Ising model with suppressed 
bulk noise [13], which we believe to be an interesting 
paradigmatic case. On the one hand as argued in Q it 
shows generalized voter like behavior, and it includes the 
voter model at its critical point. On the other hand we 
will discuss an effect on the level of single agents com- 
parable to the Sznajd model, in the presence of a social 
temperature. 



The model is defined in section [III Its dynamics is 
described briefly with special emphasis on the role of 
the model parameter as social temperature. Finally, the 
presence of absorbing states is discussed. In section IIIII 
we slightly modify the model to avoid absorbing states. 
With this stationary state model variant we can describe 
the phase transition with standard methods [14]. We find 
an abrupt jump of the order parameter, but the fourth 
order cumulant and critical divergences at the phase 
transition emphasize the continuous type of the phase 
transition. The stationary state model variant shares 
the same critical exponents as found with dynamic non- 
equilibrium methods For the stationary state model 
variant we finally establish a Fokker-Planck description in 
section HVl We compare numerical results with analytical 
mean field results and thereby find a simple explanation 
for the behavior at the phase transition. This behavior 
is also found for a small world variant of the model, be- 
ing closer to real social systems. With the Fokker-Planck 
description we also capture finite size effects. We find 
a divergence related to the absorbing states even in the 
stationary state model variant. This results in a slowing 
down of the dynamics. In section |V] we summarize our 
results and give a brief outlook. 
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II. MODEL DESCRIPTION AND SOCIAL 
TEMPERATURE 



We investigate a kinetic Ising model [13] which consists 
of N = agents on a two dimensional torus with opin- 
ions Si — ±1. Every spin adapts in random sequential 
update to his four nearest neighbors, and thus its state in 
the context of its neighborhood is adequately described 
by with the number of agreeing neighbors 



to the transition probabilities 
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jGnn(z) 



(1) 



with nn(z) denoting the nearest neighbors of agent i. Be- 
cause we want to model equivalent opinions, the proba- 
bilities for spin flips (^Si^Ui) (—5^,4 — ix^) are spin inde- 
pendent: P(n,)^(4-n,)- With the convention p(^)^(4_^) + 
P(4_^)^(^^) = 1, as in a heat bath Monte Carlo simulation 
of the Ising model, we immediately get P(2)^(2) = 1/2, 
leaving only two independent parameters p(4)^(o) (isola- 
tion, bulk noise) and p(3)^(i) (join minority, interfacial 
noise). 

Heat bath flipping probabilities for the Ising model 
(with a coupling constant of one) in our notation read 
as p(^)^(4_^) = (1 + exp(— /3(8 — 4ii))~-^ with the inverse 
temperature /3. This relates the two independent param- 
eters to each other and leads to an equilibrium system 
with defined temperature obeying detailed balance for 
single spin flips [l3| (by construction of the heat bath 
algorithm). In kinetic Ising models this relation is not 
fulfilled. Such models can be understood as being cou- 
pled to two heat baths of different temperatures [13] and 
thus they experience a flux of heat, leading to a first sign 
of non-equilibrium behavior, which will be supplemented 
by a more drastic non-equilibrium property later on. 



T=0 T<T^ T=Tv T>Tv 

|s. "^'O "^Q 



FIG. 1: Illustration of the persuasiveness an agent si receives 
from its neighbors for different temperatures T . For low tem- 
peratures, majorities convince stronger. For high tempera- 
tures, the agent follows any present opinion in a panic-like 
mood. 



P(4)^(0) = (2) 

^(0)^(4) = 1 

(no isolation / no bulk-noise) 

= l + exp(4/3) 

P(l)^(3) = 1 -P(3)^(l) 

(interfacial- noise) . 

This single parameter model for any /3 shares the prop- 
erty of the voter model, where agents only adopt opin- 
ions which actually exist in their neighborhood. This 
property is reasonable since total isolation seems to be 
a quite rare event in social systems. In the voter model 
every neighbor influences an agent with the same per- 
suasiveness, since one of his neighbors is randomly cho- 
sen for his update. Thus we get p(^)^(4_^) = (4 — 'u)/4, 
which is a special case in our model defined by p(3)^(i) = 
1/(1 + exp(4/Tv)) = 1/4 with the special temperature 
Tv = 4/ln(3) ~ 3.641. For smaller temperatures the 
model shows an increased persuasiveness of groups of 
agents, 'P{x)^{z) > ^P{3)^{i)^ as illustrated in Fig.[Tl This 
motivates the name Group- Voter Model (GRVM). This 
effect, which motivated models like the Sznajd model 
and majority models [l|, is incorporated in our model 
in a continuously tunable way. For higher temperatures 
T > Tv, local majorities have a suppressed persuasive 
power. Agents adopt any opinion in their neighborhood 
without trusting majorities which can be seen as a panic- 
like behavior. 



T<Tv 



T = Tv 




We here introduce a system without voluntary isola- 
tion of agents which means the lack of bulk-noise or in 
other words zero temperature of bulk noise. This leads 



T>Tv 



FIG. 2: Lattice states (L = 500) for T = 0.9 Ty (top), T = 
Ty (middle) and T = 1.1 Ty (bottom) at increasing times 
t = 1000 (left), t = 10000 (middle), and t = 50000 (right). 
The changing ordering dynamics indicates a phase transition 
at Ty, where the dynamics is identical to the voter model. 



In Fig. [21 the coarsening dynamics starting from ran- 
dom initial conditions can be seen for the different tem- 
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perature regimes and a system of = 500^ agents. For 
low temperatures (top line), where local majorities are 
strongly preferred and thus single agents follow group 
opinions, we see ordering dynamics leading to an ordered 
phase. So the group-following tendency has an ordering 
effect compared to the pure voter model. At T = Ty we 
get the voter model with slow cluster growth and rough- 
ened surfaces. At this temperature we find a phase tran- 
sition from order to disorder (see below). At this point 
the temporal decline of the density of surfaces changes 
from power law behavior through logarithmic decline to 
saturation i. For T > Tv we see a disordered noisy 
system, so the lack of trust of single agents in local ma- 
jorities has a destabilizing effect. 

Summarizing, an increasing temperature leads to a de- 
creasing group effect (Fig.[T]), and simultaneously it leads 
to an increased dynamical temperature (Fig.[2j). This can 
be interpreted as follows: In a faster changing environ- 
ment, any opinion provided by any neighbor seems inter- 
esting, single neighbors might be right by chance. In the 
limit of infinite temperature, opinions which are present 
in the neighborhood lead to adaption with probability 
1/2 regardless of minorities or majorities. On the other 
hand, decreasing trust in local majorities indeed acceler- 
ates the dynamics, so the described effect is self consis- 
tent. The complementary micro- and macro-effects of the 
parameter T ranging from group following and ordered 
states to panic reactions and disordered states allow us 
to interpret it as a social temperature. 

For T > Ty the system remains in a state in which the 
magnetization per site 

1 ^ 



Magnetization 



Density 



m 
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(4) 



oscillates around zero. Note that the amplitude of the 
oscillations is smaller for higher temperatures and higher 
system sizes. In Fig. [3] we finally see the presence of ab- 
sorbing states for a system with L = 75 and T = 1.02-Ty. 
The system is in the disordered state and oscillates 
around zero magnetization for more than one million 
sweeps, and the according density of magnetization 

1 v-^ ^ 



p{m') 



t 



^m' ,m{t) 



(5) 



with the normalization constant c and the Kronecker 6ij 
suggests that the system reached a stationary state. Nev- 
ertheless the system finally gets trapped in an absorbing 
state. These absorbing states exist due to the absence 
of bulk noise. In finite systems they are reached for any 
temperature and system size. Thus the system start- 
ing from random initial conditions performs a transient 
nonequilibrium dynamics. 

The description of non-equilibrium systems is some- 
what complicated. There are different approaches to deal 
with the properties associated with non-equilibrium sys- 
tems. One possible way is to investigate the phase or- 
dering dynamics and recover dynamic critical exponents 



E 




FIG. 3: (Color online) Left: Magnetization over time t for 
T — 1.02 Tv^ and L — 75. Right: Magnetization density. 
Absorbing states can be observed. 



We here use a different procedure. We slightly change 
the dynamics to prevent the system from reaching the 
absorbing state. The resulting stationary state model 
variant we will describe and investigate in the next sec- 
tion. 



III. 



STATIONARY STATE MODEL VARIANT 
AND THE PHASE TRANSITION 
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FIG. 4: (Color online) Comparison of the stationary state 
model variants with the daemon rule (left) and with bulk 
noise (a = 10~^, right) with the dashed line (daemon rule for 
L = 400) for comparison. The daemon rule reproduces the 
jump of the control parameter in the thermodynamic limit. 



Similar to [15|, we use a minimal modification of the 
original model to get an estimation of the quasi-static 
properties. We let a little daemon keep the last remaining 
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spin from flipping which prevents the system from reach- 
ing the absorbing states. The modifled model reaches a 
stationary state which is used for the evaluation of the 
statistical quantities. As shown by Dickman and Vidigal 
[14]. such methods are suitable for studying the univer- 
sal behavior of transient dynamics that otherwise would 
decay into absorbing states. As we still have a system 
coupled to two heatbaths with different temperatures, 
our system exhibits a non-equilibrium stationary state 
(NESS). 

We show that the daemon rule indeed is minimal by 
comparing it to a different procedure. In the case of a non 
zero isolation probability p^^o = aps^i with little a a 
little bit of bulk-noise is introduced. The magnetization 
for different L is shown on the right hand side of Fig. |4j 
The phase transition of the system with tiny bulk noise 
is Ising-like. The smooth transition can also be observed 
for bigger L, although a = 10~^ is already considerably 
small. This suggests that the limit a in addition to 
L ^ oo must be taken to model the phase transition of 
the GRVM. Whereas the daemon rule produces a phase 
transition which displays the abrupt change in magneti- 
zation m. This is shown on the left hand side of Fig. |4j 
Note that the impact of the daemon rule decreases for 
bigger L, whereas the impact of the method with bulk 
noise remains constant. 
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The absorbing states dominate the time series of m for 
T < Ty. For T ^ Ty the system spends a flnite time 
in states where |m| is not close to 1, but states with |m| 
close to 1 are still the most frequent states and the system 
spends a typical time of approximatively 1 million sweeps 
near one of the absorbing states. If we consider the plots 
corresponding to T = 1.0046 • Ty, we see that although a 
clear maximum in the magnetization's density p can be 
found at m = 0, the system spends a considerable time 
in a state in which |m| is close to one. This effect is a sign 
of the absorbing state around |m| = 1 due to vanishing 
diffusion, as will be discussed in the next section. The 
effect increases with system size and thus facilitates flnite 
size scaling. It must be taken into account additionally 
to the critical slow down and it dramatically enhances 
the number of sweeps needed for good statistics. 
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FIG. 6: (Color online) The fourth order cumulant Ul over 
the temperature for different system sizes . The typical 
behavior for a continuous phase transition can be observed. 



To examine the critical behavior at the phase transi- 
tion we flrst observe the fourth order cumulant Ul- We 
evaluate Ul with the time series of the magnetization 
using 



FIG. 5: (Color online) Comparison of magnetization m over 
time t (left) and its density p (right) for different temperatures 
T for L = 100. From top to bottom: T = 0.9846 Ty, T = 
1.0002 Ty, T = 1.0046 Ty and T = 1.6717 Ty. The time 
spend near the absorbing states is long for T ^ Ty, giving 
rise to a slowdown of the simulation. 



Although this method eliminates the existence of ab- 
sorbing states and produces a system with suitable sta- 
tionary state, the absence of bulk-noise gives rise to an- 
other problem. If the system is in a state where the 
daemon is needed, it is likely that the system will remain 
in this state for an extended amount of time. Time series 
of m and the corresponding densities for four different T 
with 5 million sweeps and L = 100 are shown in Fig. [5l 
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Ul can be used to determine the order of transition 
A flrst order transition would have a Ul that dis- 
plays a clear minimum. This minimum gets closer to 
the system's critical parameter as L gets bigger. A sec- 
ond order transition has a /7l that starts at | and falls 
monotonously to zero as T grows. The Ul intersects at 
the critical parameter of the system. 

The calculated fourth order cumulant (Fig. [6]) is | for 
T < Ty and falls to zero for T > Ty. This behavior is in 
agreement with a second order phase transition. The sys- 
tem has therefore a discontinuous phase transition with 
critical behavior as typical for continuous phase transi- 
tions. 
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FIG. 7: (Color online) Correlation length ^ for different sys- 
tem sizes over reduced temperature e. Dashed line: Plot 
of ^ OC to compare the data to the expected critical 

exponent. 




Since there is a jump in the order parameter {(3 = 0) and 
the dimension is d = 2, the remaining exponents obey 
1 = 2. We expect u = ^ for the universality class of the 
GVM [3]. As shown in Fig. [71 an exponent of u 
the obtained data nicely. 
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Due to the phase transition with continuous properties 
a critical behavior with scaling laws is expected. We 
evaluate the susceptibility and the correlation length to 
determine the scaling exponents 7 and u. The correlation 
length £ can be measured by evaluating the structure 
factor [13, which is often used in solid states physics. 
The structure factor is obtained using S{k) = \s\k)\'^ 
where s' = F{s'} is the Fourier transformed grid s' = 
s — m. Note that this definition ensures that states with 
large magnetization lead to small values of ^. We take 
the circular average over k and get S{k). The structure 
factor should go to zero as the magnitude of k approaches 
zero or infinity and thus S{k) displays a peak at the k' 
corresponding to the average domain size. ^ = 2i\ jk' was 
calculated by using 

^ Cfc.g(fc)dfc 

The simulations were performed with 5 million sweeps. 
After the system reached its stationary state, ^ was calcu- 
lated every 10, 000 sweeps to ensure uncorrelated results. 
An average of the domain size ^ and the error given by 
the standard deviation was therefore possible to acquire. 
The error near the critical point is relatively large. One 
could assume that the errors would be smaller if ^ was cal- 
culated more often, for example every 100 sweeps. This 
is not the case, since the discussed slowdown has to be 
considered. As shown in Fig. O the time required to go 
through a sufficient number of states is large for T ^ Ty. 
So only larger time series can improve the statistics. 

If we plot the correlation length over the reduced tem- 
perature {e = ) , we should observe the typical crit- 
ical behavior ^ ex \€\~^ . The exponents obey the relation 

j = d-u-2p. (8) 



FIG. 8: (Color online) Left: susceptibility x as a function 
of T — Ty. Right: reduced susceptibility x cis a function of 
(T - Tv)L^ = sTv with J = 2 and Tc = 0.9997 



To evaluate the remaining critical exponent 7 the sus- 
ceptibility XL is needed. We calculate xl as proposed in 
|18j by using the time series of m, 

XL=N{{m^)-{\m\f)^\er. (9) 

The exponent 7 can be found by introducing the reduced 
susceptibility given by 

x{e)=x{£Li)=XL{s)L-i. (10) 

The reduced susceptibility x plotted against e = sL^ for 
different L should collapse for T < Ty given the correct 
7, u and critical temperature T^. This was realized for 
different L (Fig. [8j). The chosen parameters to obtain 
Fig. [8] were ^ = 2 and Tc = 0.9997 Ty. 



IV. FOKKER-PLANCK EQUATION AND 
MEAN FIELD DESCRIPTION 

As a starting point for the Fokker-Planck description 
of the dynamics, we use transition probabilities for the 
magnetization p^{m) = Pm±Am,m for single time steps 
At = 1/N with Am = which describes a single pos- 
sibly flipping spin. With these two switching probabili- 
ties at hand, the Fokker-Planck equation can be deduced 
in a very simple way as described in the following. We an- 
alyze the time evolution of the probability density of the 
magnetization P(m, t) by performing a Taylor expansion 
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(AP = P(m, t + At) - P(m, t)) 
AP p'^{m — Am) 



At 



At 

p~{m + Am) 
At 

p'^{m) p~{m) 



At 
9 / Am 



P(m — Am, t) 
P(m + Am,t) 
P(m,t) 



At 



9m V At 



(j9+(m) -p-{m))P{m,t) I (11) 



1 



Am^ 



(p^(m)+p (m))P{m^t) 



2 dm? V At 
This leads to the Fokker-Planck equation 



dP 

~dt 



_d_ 

dm 



(ai(m)P(m, t)) 



1 



2dm^ 



^2 



Am^ 



(p^(m) + p (m)) 



(a2(m)P(m, t)) 
(12) 

(13) 
(14) 



with the drift term ai(m) and the diffusion term a2(m). 
Using a potential for the drift term ai(m) = — ^^(^) 
the phase transition can be understood. Additionally 
with the diffusion term a2 (m) and the stationary solution 



Ps{m) 



■ exp I 2 

«2(m) V ^mo 



ai(m') 
a2{m') 



dm' 



(15) 



with the normalization constant c, finite size effects may 
be investigated. 

To calculate the properties of interest, the transition 
probabilities p^ (m) are needed. For simulations on the 
grid, they are calculated using the time series m(t) (with 
single step resolution At = 1/A^). With Am(t) = m(t + 
At) — m(t) we get 



p'^(m) 



t <^Am(t),d=A 



m'-'m{t) ,m 



t ^m{t),m 



(16) 



We want to compare our results to mean field calcula- 
tions, which we calculate using the flipping probabilities 

P(n)^(4-n) and 



p^{m) 



u=0 



)(m) 



N 



-P(n)^(4-n) 



(17) 



with A^(s,n) being the total number of spins with the state 
(s,ii). Assuming A/'(s,n) to have a strict dependence on 
the magnetization is reasonable for the mean field case, 
where we get (including the unlikely case of twice choos- 
ing the same neighbor) 



N 

,m - 



N 



i-s) 

N 



(18) 
(19) 



and thus finally calculate 



ar (m) = -2 p(3)^(i) - T m(l - m') 



m^ 

T 



(20) 
^0 (21) 



2 

TV 



(^(1 - m^) + (^P(3)^(i) - \^ m\l - m^) 
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(22) 



exp 



1 — m^ 



r P(3)^(l) 



m^N 



(23) 



with approximate results for \pi^s)^{i) ~ 1/4| <C 1. The 
parameters (j) are needed for a fit procedure later on. 
They are (j)y = (j)a = = "^o may be chosen and c 
is a normalization constant. 



8- 10"' 
6- 10"^ 
4- 10"^ 
2- 10"^ 


-2- 10"^ 
6- 10"^ 



2- 10"^ 

1.5 
1 

0.5 




1 1 1 1 1 
- 0^=0.079 - 


1 1 1 1 1 


1 1 1 1 1 

^/ 

0^=0.094 
1 1 1 1 1 


0^=0.111 




1 III 

A 

+ 0^=0.307 + 
I \ \ — 1 


a; 

/ 0^=0.315 \ 
I \ \ — J 


1 III 

A 

/ 0,=O.331 X 
I \ \ 1 






j|. 1 1 1 4 

U 


' W=<^.29' ' 
[ ^'^^ I 

1 1 1 1 1 





-1 -0.5 0.5 1-1 -0.5 0.5 1-1 -0.5 0.5 1-1 -0.5 0.5 1 



FIG. 9: (Color online) Drift potential V{m) (top), diffusion 
term a2(m) (middle) and density p{m) (bottom) for different 
temperatures (from left to right T = 0.997 Ty, T = Tv,T = 
1.006 Tv, T = 1.016 Tv) and L ^ 100. Measured data (+) 
is compared to mean field results with good agreement. The 
fitting parameters say that diffusion and drift are suppressed 
compared to the mean field case. The drift potential clarifies 
the phase transition. The vanishing drift at \m\ = 1 leads to 
a divergence in the density, related to the absorbing states 
connected with a dynamic slowing down. 



The symbols in Fig. [9] show results for the drift po- 
tential V{m) (top), the diffusion term a2(m) (middle) 
and the density p(m) (bottom) for the following tempera- 
tures around the critical temperature (from left to right): 
T = 0.997Tv, T = Tv, T = l.OOGTv, T = l.OlGTv. 
The system size was chosen as TV = = 100^ and for 
the lower temperatures time series of 6 • 10^ sweeps were 
performed, for the higher temperatures 2-10^ sweeps suf- 
ficed. So we used altogether more than 10^^ steps in eq. 
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(p!6|) . The dashed hnes in the figure are fits of the mean 
field functions (|2T]) , (|22]) and (|23]) to the measured data 
and the fit parameters (j)a and (j)r are given (two of 
them are not defined in the critical case). As we can 
see, drift and diffusion are suppressed compared to the 
mean field case, which leads to an effect to the density, 
which might be seen as a temperature stretching proce- 
dure away from the critical temperature Ty. Apart from 
that, the system shows perfect mean field like behavior. 
With this knowledge we are able to discuss the phase 
transition including the critical exponents and the finite 
size properties of the system including the role of the 
absorbing states. 

To discuss the phase transition, we have to consider the 
case N ^ oo. From calculations performed for N = 20^ 
and N = 50^ we know, that the system behaves mean 
field like in all these cases, with shrinking parameters (/)^, 
(j)a and <pr for increasing N. The drift correction shrinks 
faster than the diffusion correction, however, diffusion is 
additionally proportional to and thus the infiuence 
of fiuctuations shrinks with increasing system size. In 
the thermodynamic limit we thus have to consider the 
minima of the potential V{m) as shown in the top line 
of Fig. [9] (see also [ll-lsl). The position of the minima 
performs a jump from \m\ = 1 to m = at T = Ty 
and thus we have an abrupt phase transition with ex- 
ponent /3 = 0. At the temperature Ty the potential 
vanishes which corresponds to the vanishing drift in the 
voter model. Adding a term —hm to the potential in eq. 
(|2T]) , we can calculate the susceptibility X = ^ by calcu- 
lating the minimum m{h) and thus prove the mean field 
exponent 7 = 1, because near the phase transition we 

get X {P(3)^(i) - ly oc (T - Tc)"^ (strictly speak- 
ing the magnetization m = ±1 in the ordered phase may 
only be infiuenced by a field h of the opposite sign). 

Finite size effects are especially visible in the densities 
as can be seen in the bottom line of Fig. [9] and in eq. (|23|) . 
respectively. Integrating the density, the term (1 — m^)~^ 
resulting from the vanishing diffusion at m = ±1 (with si- 
multaneously vanishing drift) leads to a divergence (com- 
pare [3, [HI ) . In the simulations the daemon keeps last 
spin from fiipping and thus the integration limits are re- 
duced to \m\ = 1 — 2/N. Nevertheless, the absorbing 
states can be seen around |m| = 1, as for increasing sys- 
tem size A/", the system with a temperature near the phase 
transition (and even in the disordered phase, as can be 
seen for T = 1.006 Ty in Fig. [9]) will spend an increasing 
amount of time near the absorbing state. So we again get 
a system which is trapped by the absorbing state if the 
system once moves to this state. This causes a slowing 
down in simulations, which is even more influential than 
the critical slow down in equilibrium systems and results 
in a poor convergence. 

Finally we want to discuss a small world variant of the 
system, where some of the nearest-neighbor grid connec- 
tions are rewired. It is quite meaningful to cover this case 
in the social context, since social networks are known to 
have small world properties. Because the small world 



network is somewhere in between the grid (large diame- 
ter) and the mean field case (diameter of one), and the 
results on the grid are qualitatively mean field like, the 
results for the small world variant should be mean field 
like, as well. 
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FIG. 10: (Color online) As in Fig. [9]for a small world variant 
with 5% rewired links and temperatures (from left to right) 
T = 0.998 Tv, T = Tv, T = 1.003 Ty, T = 1.006 Ty. The 
shortcuts lead to increased drift and diffusion. 



We start from the grid and use directed links, which are 
directed from infiuential agents to the infiuenced ones. So 
an agent takes into account four in-links to calculate the 
switching probabilities, or more precisely the opinions of 
the agents on the other end of the in-links. We rewire 
5% of the in-links by randomizing the starting points of 
the links and keeping the end points, which keeps all 
in-degrees to be four and only changes the out-degree 
of several agents (since single agents might have an out 
degree of zero, we have to change the daemon rule to 
act on more than one last spin). Fig. [TOl shows results 
for the small world network as in Fig. [9l here for the 
temperatures (from left to right) T = 0.998 Ty, T = Ty, 
T = 1.003 Ty and T = 1.006 Ty. We find the same mean 
field like behavior as for the grid and conclude, that the 
same kind of phase transition should appear. The fit 
parameters (j) are more close to the mean field case = 
(/>a = (/>r = 1 and we found using calculations for N = 50^ 
that these fit parameters do not remarkably change for 
increasing system size, which is in accordance to slowly 
growing diameters in small world networks. 



V. SUMMARY 

We investigate a kinetic Ising model with suppressed 
bulk noise, which is known as a prominent representa- 
tive of the generalized voter model phase transition. For 
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low temperatures, we emphasize the effect of enhanced 
persuasiveness of groups on the level of single agents. 
This leads to ordered states compared to the pure voter 
model. Similar effects are implemented in opinion for- 
mation models as for example the Sznajd model or the 
Majority Rule model. For high temperatures the oppo- 
site effect of lacking trust in majorities can be observed 
which leads to increasingly disordered states. Through 
the model parameter T the group effect can be tuned, 
which leads to a changing behavior of single agents and 
the system as a whole as well. This effect allows us to 
identify the model parameter as a social temperature. 
So we find that this single parameter system shows inter- 
esting properties not only for its behavior at the phase 
transition, but also provides intuitive rules at the level of 
single agents. 

To describe the system using its transient quasi-static 
properties, we change the dynamics using a minimal rule 
which prevents that absorbing states are reached. We 
find the generalized voter model transition, as it has been 
found for the original system [sj. With the fourth order 



cumulant we emphasize the continuous type of the phase 
transition despite the jump in the order parameter. 

Finally we derive the Fokker-Planck description of the 
phase transition. We measure drift and diffusion using 
numeric time series and compare it to analytic mean field 
results. We find perfect mean field behavior, only with 
suppressed diffusion and drift. With this we understand 
the abrupt phase transition including the value of the 
critical exponent 7. We additionally find a divergence 
in the magnetization density due to vanishing drift and 
diffusion, related to the absorbing states. This leads to 
a dynamic slowing down. 

As a system between grid and mean field, we also check 
a small world variant. This system shows exactly the 
same mean field like behavior, but diffusion and drift are 
less suppressed. This system is more similar to real social 
systems with small world properties. 

In summary we obtain a physical understanding of the 
voter model class in terms of a kinetic Ising model, obey- 
ing mean field dynamics. 



[1] C. Castellano, S. Fortunate, and V. Loreto, Rev. Mod. 

Phys. 81, 591-646 (2009). 
[2] P. Clifford and A. Sudbury, Biometrika 60, 581-588 

(1973). 

[3] I. Dornic, H. Chate, J. Chave, and H. Hinrichsen, Phys. 

Rev. Let. 87, 045701(2001). 
[4] H. Hinrichsen, Advances in Physics 49, 815-958 (2000). 
[5] G. Odor, Rev. Mod. Phys. 76, 663 (2004). 
[6] O. Al Hammal, H. Chate, L Dornic, and M.A. Munoz, 

Phys. Rev. Let. 94, 230601 (2005). 
[7] C. Castellano, M.A. Munoz, and R. Pastor- S at orr as, 

Phys. Rev. E 80, 041129 (2009). 
[8] F. Vazquez and C. Lopez, Phys. Rev. E 78, 061127 

(2008). 

[9] K. Sznajd- Weron and J. Sznajd, Int. J. Mod. Phys. Cll, 
1157-1165 (2000). 



[10] S. Galam, Eur. Phys. J. B 25, 403406 (2002). 

[11] S. Bornholdt, Int. J. Mod. Phys. C 12, 667-674 (2001). 

[12] S.M. Krause and S. Bornholdt, arXiv: 1103. 5345 vl. 

[13] M.J. de Oliveria, J.F.F. Mendes, and M.A. Santos, J. 

Phys. A 26, 2317 (1993). 
[14] R. Dickman and R. Vidigal, J. Phys. A 35, 11471166 

(2002). 

[15] M.A. Munoz, Phys. Rev. E 57, 1377-1383 (1998). 

[16] D.P. Landau and K. Binder, A guide to Monte Carlo 
Simulations in Statistical Physics, Cambridge University 
Press, New York (2000). 

[17] M.E.J. Newman and G.T. Barkema, Monte Carlo Meth- 
ods in Statistical Physics, Oxford University Press, Ox- 
ford (2004). 

[18] M.J. de Oliveira, J. Stat. Phys. 66, 273 (1992). 



